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ABSTRACT 

We present results obtained from a series of 5 A resolution spectra of 
the X-ray Nova GRO J0422-F32 obtained in 1993 October, when the system 
was approximately 2 magnitudes above quiescence, with R ~ 19. The data 
were obtained in an effort to measure the orbital radial velocity curve of the 
secondary, but detection of the narrow photospheric absorption lines needed 
to do this proved elusive. Instead we found wide absorption bands reminiscent 
of M star photospheric features. The parameters determined by fitting 
accretion disk line profiles (Smak profiles) to the Ha line are similar to those 
found in several strong black-hole candidates. Measurements of the velocity 
of the Ha line are consistent with an orbital period of 5.1 hours and a velocity 
semi-amplitude of the primary of 34 ± 6 km s~^. These measurements, when 
combined with measurements of the velocity semi-amplitude of the secondary 
made by others, indicate that the mass ratio g ~ 0.09. If the secondary 
follows the empirical mass-radius relation found for CVs, the low q implies 
a primary mass of ~ S.GM©, and a rather low (face-on) inclination. 



^Data reported herein was obtained with the Whipple Observatory MMT, an instrument jointly 
operated by the Smithsonian Astrophysical Observatory and the University of Arizona. 
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The Ha EW is found to be modulated on the orbital period with a phasing 
that implies a partial eclipse of the disk by the secondary, but simultaneous 
R band photometry shows no evidence for such an eclipse. 



1. Introduction 

Low Mass X-ray Binaries (LMXBs) are systems where a neutron star or black hole 
accretes material from a low mass (M<1M0) secondary (sec Bhattacharya and van 
den Heuvel 1991). X-ray novae (XRN) form a subset of LMXBs whose properties are 
characterized by dramatic increases in flux (factors of ~10^-10^) at X-ray and radio 
energies on time scales of days, with a subsequent decay to quiescence on time scales of 
months. 

While these transients are in outburst, their X-ray and optical properties arc very similar 
to those of the persistently bright LMXBs, where the optical flux is dominated by the 
reprocessed emission from the X-ray irradiated disk. In quiescence, however, the disk 
often fades to reveal the secondary star itself. This allows detailed photometric and 
spectroscopic measurements of the secondary to be made - otherwise impossible in the 
outburst state, and in most of the persistently bright LMXBs. 

Studies of quiescent XRN arc especially important for the following reasons. The 
frequency of black hole candidates among X-ray novae appears to be remarkably high; 
9 out of the 13 optically identified systems are strong or probable black hole candidates 
(van Paradijs and McClintock 1995). Our understanding of the evolution of LMXBs (eg, 
the progenitors, lifetimes, and numbers) will be limited until we understand the nature 
of the transient sources. This is because the number of "dormant" transients may well 
exceed the number of persistent LMXBs by one to two orders of magnitude (van den 
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Heuvel 1992). 

The detection of the photospheric hnes of the secondary in quiescent systems allows 
measurement of the orbital radial velocity curve, leading to firm lower limits on the mass 
of the accreting object. These measurements indicate masses t 3Mq for the X-ray nova 
A0620-00 (McClintock and Remillard 1986), V404 Cyg (Casares, Charles and Naylor 
1992) Nova Muscae 1991 (Remillard, McClintock and Bailyn 1992), Nova Ophiuci 1977 
(Remillard et al. 1995), GRO J1655-40 (Bailyn et al. 1995), and GS2000+25 (Casares, 
Charles and Marsh 1995). Because the maximum mass for a stable neutron star is 
~ 3Mq (Chitre and Hartle 1976) these six XRN are regarded as strong black hole 
candidates. The mass ratio q can be determined by combining the measurements of the 
orbital velocity of the secondary with similar measurements of emission lines originating 
in the accretion disk. Models for formation of emission lines in accretion disks depend 
on many parameters, among them the inclination, inner and outer disk radii, and central 
object mass, and comparison to observed profiles can help to constrain some of these 
parameters. 

In Callanan et al. 1995 (Paper I) we discussed the discovery of GRO J0422-I-32, and 
presented the first results from our optical observing campaign during the outburst and 
decline. In this paper we present the results of extensive spectroscopic observations 
obtained ~ 430 days after outburst, when GRO J0422-I-32 had decayed to R 19 (see 
Figure 1). Our continuing photometric observations show that the source subsequently 
reached its quiescent level at day ~ 760. The average intensity during days 760-950 
is R = 20.94 ± 0.11 (28 measurements) and V = 22.35 ± 0.17 (3 measurements). We 
note that the error on the R magnitude in quiescence quoted here is computed from the 
observed scatter in the data, and may be dominated by intrinsic variability in the source. 
The outburst amphtude was AV ~ 9.0, which is the largest yet seen for an XRN: typical 
outburst amplitudes are AV = 5 — 7 magnitudes (van Paradijs and McClintock 1995). 
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2. Observations 

By 1993 October 10, ~ 14 months after the initial outburst, GRO J0422+32 had reached 
R ~ 19 (day 430 on Figure |l]), but was still at least ~ 2 magnitudes brighter than its 
quiescent level. This is surprising given that the black hole candidates Nova Muscae 1991 
(Remillard, McChntock and Bailyn 1992), GS2000+25 (Chevalier and Ilovaisky 1990), 
and V404 Cyg (Casares and Charles 1992) had reached quiescence within 14 months of 
the outburst. 

Spectra of GRO J0422+32 were obtained with the SAO Whipple Observatory MMT Red 
Channel Spectrograph on 1993 October 10 and 12 (UT), and spectra of several bright 
standards were obtained on October 11 (UT) for use in determining absolute radial 
velocities and spectral types. The spectra were obtained with a Loral 1200 x 800 pixel 
CCD having 7 electrons read noise. The combination of the grating and 1" slit gave 
~ 5 A FWHM spectral resolution and spectral coverage from ~ 5000 A through 
~ 6700 A. The seeing on October 10 and 12 was predominantly ~ 1" and the skies were 
mostly clear, but occasional clouds and intervals of worse seeing (~ 2") compromised 
the photometric accuracy of our spectra. The cloud cover on October 11 was sufficiently 
heavy that only bright standards were visible on the TV guider; however GRO J0422+32 
was easily visible on the other two nights. The slit was approximately maintained at 
the parallactic angle in order to maximize throughput. Twenty minute integrations of 
GRO J0422+32 were bracketed by short exposures of a HeArNe lamp for wavelength 
calibration purposes. Data were reduced and analyzed using standard IRAF routines. 
Sixteen spectra were obtained on each night for a total exposure time of 640 minutes 
on GRO J0422+32. Simultaneous Whipple Observatory 1.2 m R band photometry was 
obtained; the instrumentation and data reduction are described in Callanan et al. 1995. 
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3. Analysis and Results 
3.1. Cross-Correlation Results 

Our primary goals were to detect the expected photospheric absorption lines due to the 
presumed late-type companion, and to measure the Doppler shifts in those lines due 
to the expected orbital motion of the star. We attempted this by cross-correlating the 
32 individual GRO J0422-I-32 spectra against the spectra of late-type dwarf templates, 
ranging from early K to mid M, which were obtained with the identical instrumental 
setup on Oct 11. Before cross-correlating, the regions around the 5577 A and 5890 A 
night-sky lines were replaced with a fit to the local continuum, because these night sky 
lines are strong enough to render subtraction in the individual 20 minute exposures 
problematic. In addition, the spectra were filtered in order to remove low-frequency 
features and in order to match the instrumental resolution of 5 A. The low frequency 
filtering was achieved by subtraction of a 13*^ order cubic spline fit to the continuum, 
and removal of wavenumbers < 6 from the Fourier transform of the spectra (FT). This 
effectively removes all features on scales larger than ~ lOOA, including the instrumental 
response and possible TiO bands (see Section 3.2). The high frequency filter removed 
wavenumbers > 330 from the FT, which has the effect of smoothing the data on the 
intrinsic resolution of 5A (~ 3.5 pixels per resolution element). In order to avoid the 
contaminating effects of the H-Balmer emission lines, only the region between 5000 A 
and 6500 A was used for the correlations. 

We were not able to detect any narrow, late-type photospheric lines in this data 
set, despite extensive efforts. In order to quantify our lack of detection, we added 
various amounts of a late-type stellar spectrum (with NaD excised) to our individual 
GRO J0422-I-32 spectra, and computed cross-correlations. In keeping with the short 
period, we added dwarf (rather than giant) spectra. When a M3V spectrum accounts 
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for ~ 10%, a MOV for ~ 20%, or a K3V for ~ 30% of the total light, the correlations 
yield significant detections with R (Tonry and Davis 1979) values in excess of 3.0. In 
this way we estimate that the fractional contribution of any late-type stellar photosphere 
to the continuum light at R ~ 19 is less than 10% ± 5%, 20% ± 5%, or 30% ± 5% for an 
assumed M3V, MOV or K3V companion, respectively. 



3.2. Shape of the Continuum 



Figure ^ shows the sum of the spectra obtained on October 10 (bottom) and 
October 12 (top). The most obvious features are double peaked Ha emission and 
interstellar NaD 5890 A absorption. The NaD night sky line has not been excised 
because it is accurately subtracted in the summed spectrum, while the 5577 A has once 
again been removed. The NaD EW (1.2 A) is consistent with an interstellar origin. Less 
obvious in Figure ^ are troughs in the continuum with minima at ~ 6200 A, ~ 5900 A, 
and ~ 5200 A. These troughs are more apparent in Figure |^, which shows a scaled 
plot of the October 10 GRO J0422+32 spectrum and a spectrum of a M3V star. The 
correspondence in the broad continuum features leads us to suspect that something like 
an M star photosphere is present in the spectra, even though cross-correlations show 
that the narrow lines expected from an M star are not detectable. We considered the 
possibility that these troughs might be an artifact of our data acquisition or reduction 
procedure; however, spectra of the nova V849 Oph {R ~ 18) which were taken on the 
same nights with the identical instrumental setup and reduced in the same way do not 
show these features. In addition, spectra taken by other observers a few months later also 
show evidence for bands indicative of a late-type stellar photosphere (Bonnet-Bidaud 
and Mouchet 1994), although the lower S/N and resolution of these spectra make the 
bands difficult to detect. The evidence for these bands was much weaker in our spectra 
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from October 12. This may be because the S/N in the October 12 data was shghtly 
lower than on October 10 (see Figure |]), or it may be that the band strength is variable 
and the ~ 0.3 magnitude brightening on October 12 is related to their weaknesses. We 
therefore concentrate solely on the October 10 spectra for the remainder of this section. 

We measured the EW of the strongest two features (at 5900A and 6200A) and compared 
them to the EW in the M3V star spectrum. We choose to compare to an M3V 
star because it is a good match to the spectral type of GRO J0422+32 determined 
in quiescence, which is M2±1V (Garcia et al. 1994, Filippenko et al. 1995). The 
continuum points and band limits are defined by the M3V star, and are 5805A-6039A 
and 6145A-6350A, as marked on Figure 3. The EW of these bands in GRO J0422+32 
are ~ lOA, while they are ~ 60A in the M3V star. Thus we see that 1/6, or ~ 17% of 
the hght may be due to an M3V star. 

An independent estimate of the amount of apparent M star photosphere present in the 
spectra can be made by creating synthetic spectra consisting of the sum of an M3V star 
and an accretion disk, and looking for the best match to the GRO J0422+32 spectra. 
This method is sensitive to the slope of the continuum, whereas the previous method 
is sensitive only to the depth of the bands. We first de-reddened the GRO J0422+32 
spectra for Ay = 1.2 (see Callanan et al. 1995), and then sum variable amounts of an 
M3V star and an F\ = F * X°' accretion disk spectrum. We choose a power law for the 
accretion disk spectrum because it is a sufficiently accurate representation of what has 
been seen in quiescent XRN. We measured the slope of the powerlaw from the published 
data on three XRN in quiescence. The results are shown in Table 1, along with the slope 
measured for GRO J0422+32 during the decay by Shrader et al. 1994. We then assumed 
the slope for GRO J0422+32 was equal to the mean of the values in Table 1, but in 
order to determine the errors below we allowed a to vary over the full range shown in 
the table. 
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We then searched for the minimum between the observed and synthetic spectra w.r.t. 
the fraction (F) of accretion disk. The minimum was not found to be unity, in 
part because the observed spectrum does not accurately match the narrow hues of the 
template M star. In order to estimate an error range for our fit, we re-scaled the error 
bars so that the best fit = 1 (see Shahbaz et al. 1994), and computed 68% confidence 
intervals at xLin + (Avni 1976). In the best fit synthetic spectrum, the M3V star 
contributes between 20lj;4% 

5OOOA and 50^28% at 6500A, where the error ranges 
include the a contribution due to the range in slopes found in Table 1. 

Given the rather substantial errors associated with measuring the contribution based 
on the slope of the continuum, we see that both measurements are consistent with 
the assumption that an M3V star photosphere contributed ^20% of the light in 
October 1993, when R ~ 19. This is marginally consistent with our lack of a result from 
the cross-correlations, as we would need ~ 20% of the light to be due to a late-type 
stellar photosphere in order to obtain reliable cross-correlation results. 



3.3. Ha profiles 

The double peaked Ha profiles shown in Figure |^ are similar to those seen in cataclysmic 
variables, where they are believed to arise because of Doppler shifts in an accretion disk 
around the compact object (Smak 1981, Home & Marsh 1986). In order to test this 
hypothesis for GRO J0422+32 we have fit model profiles appropriate for optically thin, 
flat, Keplerian disks. The model and fitting methods are described in detail by Smak 
1981 and Orosz et al. 1994. We restricted ourselves to models in which a = 1.5 (a is the 
exponent of the radially dependent emission function), as this value is generally seen in 
dwarf nova disks in quiescence (Horne 1993), and has also been found in XRN disks in 
quiescence (Johnson, Kulkarni and Oke 1989, Orosz et al. 1994). If we use other values 
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for a the results of our model fits might be different (see Figure 4 of Johnson, Kulkarni 
and Oke 1989), but by restricting ourselves to a = 1.5 we allow direct comparison to 
previous work, in particular Orosz et ai. 1994. 

As the steady state Smak model which we fit to the data predicts a symmetrical emission 
line profile, we begin by selecting 17 of our 32 spectra from October 10 and 12 such that 
the sum has red and blue peaks of equal amplitude. A sum of all 32 spectra produces 
a profile in which the blue peak dominates (Figure This departure from symmetry 
might be explained by phase or time variable emission from an accretion disk hot spot 
seen in CVs and XRN (Johnson et ai. 1989), or it might be due to poorly understood 
deviations from optically thin, Keplerian models which have previously been noted in 
XRN (Orosz et al. 1994, Haswell and Shaffer 1990). We have not attempted to correct 
the average Ha profile for the possible effects due to the motion of the compact primary 
or the Ha absorption line expected from the secondary, as both of these effects are 
expected to be small in GRO J0422+32 and are typically ignored (ie, Johnson, Kulkarni 
and Oke 1989, Orosz et ai. 1994). In any case, fits to the simple, symmetrical Smak 
models will allow comparison of the GRO J0422+32 line profile to those seen in other 
XRN and CVs (Orosz et ai. 1994, Williams 1983.) 

A first-order polynomial was fit to the 6350 A - 6750 A region (excluding Ha) in order 
to subtract the continuum from the summed spectrum, and a chi-squared minimizing 
method was used to determine the best fit model parameters. The best fit parameters 
and 68% confidence intervals that we found are: ri = 0.14 ± 0.01, Ao = 6564.0 ± 0.1 A, 
and Vd = 442 ± 10 km s^^, were ri is the ratio of the radii of the inner and outer edges of 
the Ha emitting part of the disk, Aq is the central wavelength of the profile, and is the 
velocity at the outer edge of the disk (which is very nearly equal to 1/2 the separation of 
the peaks, see Smak 1981). A comparison of the data and best fit profile can be seen in 
Figure ^. 



3.4. Ha Radial Velocity 



Because the Ha line forms in the accretion disk, one might expect that the 
line velocity reflects the motion of the compact object. Measurement of 
this velocity should then determine the orbital period and mass function, 
/(m) = PKI/2t:G = (McSin(i))7(M^ + Mcf, where P is the orbital period, 
is the semi- amplitude of the orbital velocity variation of the primary (and of the Ha 
line by our assumption), Mc is the mass of the companion, the mass of the compact 
star, and i is the orbital inclination. The phase and mid-point (7 velocity) of the Ha 
velocity curve relative to that of the secondary provide checks on the validity of the 
assumption. Studies of CVs (eg. Young, Schneider and Shectman 1981, Stover 1981) and 
XRN (Haswell and Shaffer 1990 hereafter HS90, Orosz et al. 1994, Marsh, Robinson, 
and Wood 1994, hereafter MRW94) show that determined from Ha emission lines is 
indicative of the velocity of the primary, even though the phase and 7 of the Ha velocity 
curve are sometimes slightly offset from that expected. 

However, the complex and variable shape of the Ha line makes it difficult to measure a 
velocity that reliably traces the motion of the compact object. Several methods have 
been developed to help alleviate this problem, among them "double Gaussian fitting" 
(Schneider and Young 1980, Shafter 1983, Marsh ct al. 1994) and Smak profile fitting 
(Orosz et al. 1994). We have employed a variant of these methods. A single Gaussian is 
fit to the outer (red and blue) portion of the Ha line, excluding the central part of the 
line which can be affected by bright spot emission, optical depth effects, and other effects 
which may not reflect the motion of the compact object. In practice this outer portion is 
determined to be that portion in which the second derivative of the flux vs wavelength 
is > as one approaches the hue center. Because this is the section of the line in which 
the profile is steepest (the derivative is highest), its location is very sensitive to changes 
in the line velocity. We tested our method on data obtained for XRN Muscae 1991 by 
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Orosz and collaborators in February 1993, and we find values of Kr^ and Tq (but not 7) 
that are consistent with those derived by the double Gaussian fitting and phase resolved 
Smak profile fitting methods (Orosz et al. 1994). 

In order to improve the S/N and search for radial velocity variations the individual 
20 minute GRO J0422+32 spectra were grouped into 6 equal phase bins (with 4-7 spectra 
in each bin) based on the suggested period of P = 5.0944 ± 0.0017 hours, (Chevalier and 
Ilovaisky 1994b; we chose this period because it has the smallest error among the many 
consistent periods reported to date). For each profile a minimizing method was used 
to determine the central wavelength, width, and amplitude of a fitted Gaussian and linear 
continuum model. The errors in the fitted parameters were determined by the rms to 
the fit. A sine curve with period fixed at 5.0944 hours was fitted to the velocities, which 
gave an amplitude = 34 ± 6 km s~^and 7 = 142 ± 4 (see Figure ^). For comparison, 
Filippenko, Matheson and Ho (1995) find = 42 ± 2 km s^^and 7 = 26.0 ±1.6 from 
recent observations in quiescence with the Keck telescope. We adopt the phase convention 
of Orosz et al. 1994 and Johnson et al. 1989, in which Tq marks the closest approach of 
the compact object to the observer, and find Tq = HJD 2, 449, 270.978 ± 0.004. 



3.5. Ha EW and Continuum Variations 

Previous studies of interacting binaries have shown that the relative phasing of variations 
in emission line velocities, emission line EW, and wide-band photometry has the potential 
to reveal much about the geometry of a binary system (eg, Pringle and Wade 1985, 
Mauche 1990). In order to measure the EW of the Ha emission line, we fitted a low-order 
polynomial to the region around the emission line in each of the 32 spectra, in order to 
determine the continuum level. Figure |^B shows the resulting Ha EW vs the phase we 
determined above. We find that the EW is modulated by a factor of ~ 2 on the orbital 
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period, varying between 20 A and 40 A. There is a broad minimum in the EW around 
phase 0.5. 

We obtained 3 nights of R band photometry with the SAO 1.2 m telescope during 
October 1993. Two of these nights were simultaneous with the MMT spectra 
described herein, and are plotted vs phase in Figure ^p. These data are also shown 
in Callanan et al. 1995, Figure 3(d) [but note that they are labeled via MST dates, 
ie, October 9 and October 11]. The first seven measurements on October 12 are 
significantly (~ 0.2 magnitudes) higher than the rest, indicating a flare-like event from 
GRO J0422+32. We excluded these points from the analysis below. The solid line is the 
average magnitude in each of 12 phase bins. This average varies by 12%, and the mean 
rms in each of the 12 bins is 5%. There appears to be a broad maximum between phases 
0.5 and 1.0, with a slight indication of a local maximum around phase 0.5. 



4. Discussion 



The apparently periodic modulation of the Ha EW and velocity reported herein agree 
with the suggested ~ 5.1 hour orbital period for GRO J0422+32 (Chevalier and Ilovaisky 
1994b), which has recently been confirmed by spectroscopic observations in quiescence 
(Orosz and Bailyn 1994, 1995; Filippenko et al. 1995, Casares et al. 1995). 

The lack of a detection of the secondary from our cross-correlation efforts indicates 
that the secondary contributed only a small fraction of the light during October 1993. 
The fact that the system faded another 2 magnitudes since October 1993 by itself 
sets an upper limit of 16% to the contribution of the secondary in October 1993. The 
non-detection is consistent with tests of the sensitivity of the cross-correlations to 
detecting the secondary in the individual 20 minute spectra, which show that we would 
need a contribution of ~ 20% for reliable detection (Section 3.2). Somewhat surprisingly 
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the secondary does seem to be discernable in the grand sum spectrum, in that wide 
(but weak) M-star bands are seen in the spectrum. The measured strength of the bands 
may be somewhat larger than the 16% allowed contribution, however the errors in the 
measurement are large, and given the tendency to overestimate the strength of weak 
absorption features in noisy spectra it seems very likely that the bands are due to the 
M2V secondary. 

The accretion disk parameters (and la errors) derived using the Smak profile fits for 
GRO J0422+32 are compared to those derived for the XRN A0620-00 and X-ray Nova 
Muscae 1991 (Orosz et al. 1994) in Table 2. We see that the ri values derived for 
A0620-00 and GRO J0422+32 are very similar, and are a factor of two larger than that 
found for Nova Muscae. Interpreted in the context of the Smak model, this indicates 
that the ratio of the inner to outer radii of the Ha emitting part of the disk is ~ 7 in 
A0620-00 and GRO J0422+32, but ~ 14 in Nova Muscae. 

We can see in Figure | that the Ha EW decreases by ~ 50% at phase 0.5, about the 
same time that the R band photometry increases by ~ 12%. The relative amplitude 
of the EW and R band variations shows that only a small part of the the Ha EW 
variations can be due to variations in the continuum. Casares et al. (1995) also find a 
~ 50% orbital modulation in the Ha EW of GRO J0422+32 in data taken in quiescence. 
This variation is in marked contrast to that seen in A0620-00 (MRW94) and V404 Cyg 
(Casares and Charles 1992, Casares et al. 1993), in which the Ha EW modulation is 
largely explained by modulation in the underlying continuum. The Ha velocities imply 
that the secondary is closest to the observer when the EW is smallest, which suggests 
that the EW variations may be due to an eclipse of the disk by the secondary. However, 
this is untenable for a simple azimuthally symmetric disk. The problem is that any 
eclipse would have to cover ~ 50% of the Ha emitting region while at the same time 
the overall R band continuum (which must be due almost entirely to the disk) peaks. 
It is perhaps more likely that the EW variations are due to variations in the projected 
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geometry and integrated emissivity of the disk itself, for example due to an azimuthally 
varying disk thickness. 

We searched for, and found, orbital radial velocity variations in the Ha emission line. 
Such variations often measure the motion of the compact object in SXT (Orosz ct al. 
1994, HS90, MRW94). In the particular case of GRO J0422+32, we can check to see if 
this velocity variation is consistent with the compact object by comparing the 7 velocity 
and phase to those measured for the secondary. 

We find a central velocity of 7 = 142 ± 4 km s"^ for the Ha hue in GRO J0422+32, 
which is not consistent with velocity of the secondary of 7 = 10 ± 4 km s~^(Filippenko 
et al. 1995). However, the 7 velocities derived from SXT Ha studies are somewhat 
problematic. The values derived depend strongly on the details of the function used to 
fit the hue (ie, HS90, Orosz et al. 1994). It is therefore not surprising that different 
methods lead to different results. In the best studies system, A0620-00, Orosz et al. 
find 7 = 78 ± 3 via Smak profile fitting, but 7 = 1.5 ± 0.8 via double Gaussian fitting. 
Even when consistent methods are used, the Ha 7 in A0620-00 appears to change from 
year to year, perhaps because of variations in the shape of the Ha profile (using double 
Gaussian fitting, HS90 find 7 = 28 ± 6 in Dec 89/Jan 90, Orosz et al. find 7 = 1.5 ± 0.8 
in Jan 91/Feb 91, and MRW94 find 7 = ± 5 in Dec 91/ Jan 92). One is lead to the 
conclusion that Ha 7 velocities are not robust measures of the binary systemic velocity, 
and that in the case of GRO J0422+32 our measurement of 7 = 142 ± 4 km s~^ argues 
neither strongly for nor against the assumption that the Ha velocity variations (K^) 
are representative of the compact object. However, even given the problematic nature 
of the Ha 7 velocities, it appears that Ha velocities are accurate, in that they yield 
results which agree with those determined by independent methods (Orosz et al. 1994, 
MRW94). 

The orbital period of GRO J0422+32 is not presently known with sufficient accuracy 
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to compare the phase of the secondary radial velocity curve measured in quiescence 
(Filippenko et al. 1995, Casares et al. 1995) to the Ha curve we measured in October 
1993. However, the large Ha EW variations, seen once per orbital cycle, might provide 
an accurate fiducial mark. If one assumes this variation is stable with respect to the 
orbital phase, then we can use the EW curve to compare our Ha velocities to those 
measured for the secondary. Comparing our Figure 5 A/B to Figure 1 of Casares et al. 
(1995), we see that the peak of the Ha EW modulation occurs when the Ha velocity 
crosses from blue to red, and the secondary velocity crosses from red to blue. This is the 
phasing expected if the Ha velocities track the compact object. 

By combining our measurement of the velocity semi-amplitude of the primary 
{Kx = 34 ± 6 km s~^) with the recently determined velocity semi- amplitude of the 
secondary {K^ = 380 ± 6 km s~^, Filippenko et al. 1995) we can determine the 
mass ratio q = K^/Kc = 0.089 ± 0.016. This mass ratio, while extreme, is similar to 
that found in other XRN and reinforces the idea that XRN as a class may have mass 
ratios extreme enough to excite SU UMa type superhumps in outburst (Callanan and 
Charles 1991, Bailyn 1992). Indeed, superhumps during outburst have been reported for 
GRO J0422+32 (Kato, Mineshige and Hirata 1995). 

The mass function sets a lower limit to the mass of the compact object of 
Mx — 1.21 ± O.O4M0 (Filippenko et al. 1995). When we use our estimate of the 
mass ratio (and the definition of the mass function) we find the lower limit is raised 
to Mj. > 1.43 ± 0.09 sin~^{i)MQ. The corresponding lower limit to the mass of the 
secondary is > 0.12 ± 0.02 sin-^{i)MQ. 

A reasonable upper limit on the mass of the secondary might be derived from the 
empirical mass/period relation found for secondaries in CVs by Patterson (1984). We feel 
this is an upper limit because studies of Cen X-4, for example, show that the secondary 
in this X-ray binary is substantially undermassive (Remillard and McClintock 1990). For 
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the 5.1 hour period of GRO J0422+32, Patterson's relation predicts = O.SMq. Given 
the q found above, this translates into an upper limit for the mass of the compact object 
of Mj, < 5.6 ± 1.4Mq, and (through the mass function) a lower limit to the inclination 
of 35°. The lack of optical and x-ray eclipses sets an upper limit to the inclination of 
i < 80° at the observed mass ratio (Chanan, Middleditch, and Nelson 1976). The rather 
modest orbital modulation of the lightcurve compared to that seen in other XRN favors 
the lower end of the allowed range in i, and therefore the higher end of the allowed range 
in Mj. We note that the q and we find are consistent with those determined from 
the observations of super humps during outburst by Kato et al. (1995). However detailed 
modehng of the hghtcurve (eg Haswell et al. 1993) will be necessary to further constrain 
i, q, and M^. 
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Table 1: 



XRN 


Wavelength Range 


a 


Reference 


A0620-00 


3000A- lOOOOA 


-1.76 


Oke 1977 


XN Muscae 1991 


5000A- 6500A 


-1.88 


Orosz et al. 1995 


V404 Cyg 


4400A- 7000A 


-1.63 


Casares et al. 1993 


GRO J0422+32 


4000A- 7000A 


-1.80 


Shrader et al. 1994 



Table 2: 



parameter 


A0620-00 Nova Muscae J0422+32 


ri 

Vd 


0.15 ±0.01 0.07 ±0.01 0.14 ±0.01 
550 ± 10 450 ± 10 442 ± 10 
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FIGURE CAPTIONS 

Figure 1: The long term light curve of GRO J0422+32 in the R band, beginning a 

few days after discovery in 1992 August and continuing for more than two years. The 
spectroscopic observations reported here took place during days 430 and 432, in between 
the two mini-outbursts (at days ~ 370 and ~ 500) which where were detected during 
the decay. The relative constancy of the data after day ~ 760 indicates that the system 
has now reached its quiescent magnitude of R = 20.94 ±0.11. An earlier version of this 
lightcurve appeared in Callanan et al. 1995. 

Figure 2: The grand sum of the spectra obtained on 1993 October 10 (bottom) and 
October 12 (top). The y-axis is flux in units of ergs cm~^sec~^A~^, with an offset 
of 0.2 added to the October 12 spectrum in order to clarify the plot. Simultaneous 
photometry from the 1.2 m telescope indicated that the mean brightness was slightly 
lower on October 10 and higher on October 12, with R — 19.2 and R — 18.9, respectively. 
The most prominent feature in the spectra is the double peaked Ha emission line. These 
spectra have not been de-redenned. The EW of the NaD line (1.2 A) is consistent with 
an interstellar origin. 

Figure 3: Summed spectrum of GRO J0422+32 from 1993 October 10 compared to a 
M3V spectrum obtained as a cross-correlation template. The GRO J0422-I-32 spectrum 
has been de- reddened by — 1.2. The straight lines indicate continuum points used to 
measure EW of absorption bands. 

Figure 4: The summed Ha profile of GRO J0422+32 and best fitting Smak model. The 
Smak model is appropriate for optically thin, Keplerian accretion disks. 
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Figure 5: Panel A shows the wavelength of the center of the Ha emission line as 
observed at the telescope and as determined from Gaussian fits to the steepest section of 
the line, assuming a 5.0944 hour orbital period, plotted vs phase. The best fitting sine 
curve of = 34 ± 6 km s~^and heliocentric 7 = 142 ± 4 km s~^ is superposed. The 
time of closest approach of the primary is Tq = HJD2, 449, 270.978 ± 0.004. Panel B 
shows the EW of Ha vs the phase we have determined from the Ha velocities. Panel C 
shows the R magnitudes from the SAO 1.2 m telescope, measured simultaneously with 
the Ha EW. The circles are from October 10, the triangles from October 12, and the 
crosses are the "flare" data (see text) from the beginning of October 12 offset by 0.2 
magnitudes. The solid line is the mean flux in each of 12 phase bins, and single error bar 
of 0.05 mag, representing the mean rms of the binned data, is shown. 
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